Trapping reactions with subdiffusive traps and particles characterized by different 

anomalous diffusion exponents 
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A number of results for reactions involving subdiffusive species all with the same anomalous 
exponent 7 have recently appeared in the literature and can often be understood in terms of a 
subordination principle whereby time t in ordinary diffusion is replaced by f . However, very 
few results are known for reactions involving different species characterized by different anomalous 
diffusion exponents. Here we study the reaction dynamics of a (sub) diffusive particle surrounded by 
a sea of (sub) diffusive traps in one dimension. We find rigorous results for the asymptotic survival 
probability of the particle in most cases, with the exception of the case of a particle that diffuses 
normally while the anomalous diffusion exponent of the traps is smaller than 2/3. 
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INTRODUCTION 



In the traditional version of the trapping problem, a 
normal diffusive (Brownian) particle (A) wanders in a 
medium doped at random with static traps (B), and dis- 
appears when they meet, A+B — > B. This problem dates 
back to Smoluchowski's theory of reaction rates at the be- 
ginning of last century, and is one of the most widely in- 
vestigated and applied problems of non-equilibrium sta- 
tistical mechanics QjEMIS- An important variation of 
the basic trapping problem, in which a diffusive particle 
wanders in a medium in which the traps are also diffu- 
sive, has been the subject of intense research since the 
seminal work of Toussaint and Wilczek . 

The principal quantity of interest in the trapping prob- 
lem is the survival probability P{t) of the A parti- 
cles. From this survival probability one is able to cal- 
culate essentially all other quantities of practical inter- 
est. Yet this probability is usually difficult to calcu- 
late, and the few instances in which it has been obtained 
are considered landmark contributions. In 1988, Bram- 
son and Lebowitz 6, 7] proved rigorously that the long- 
time survival probability of a particle diffusing in a one- 
dimensional medium doped with diffusive traps decays 
as P{t) ~ exp(— At^/^), A being an undetermined param- 
eter. The evaluation of this constant proved elusive for 
many years, engendering much confusion and proposed 
solutions that were mutually contradictory. Finally, quite 
recently Bray and Blythe jSj proved in a simple and 
elegant way, assuming the so-called "Pascal principle" 
(see below), that the survival probability of a diffusing 
particle with diffusion coefficient D' in a d-dimensional 
medium with d < 2 in which the traps are also diffusive 
with diffusion coefficient D is independent of D' for long 
times. They furthermore proved that the survival proba- 
bility coincides with that of an immobile target {D' = 0) 
in the presence of diffusive traps. In particular, in a one- 
dimensional medium P{t) ~ exp(— 4p(Z?i/7r)^/^), where 



p is the density of traps. Bray and Blythe obtained their 
results by calculating an upper and a lower bound for the 
survival probability that converge to one another asymp- 
totically. Some (but not all) of the bounding results of 
Bray and Blythe have been extended by Oshanin et al. (9j 
to systems where the traps perform a compact explo- 
ration of the space, i. e., where the fractal dimension 
dui of the trajectories of the traps is greater than the 
dimension d of the space. 

The Pascal principle states that the best strategy for 
survival is for the A particle not to move. This assump- 
tion was adopted in one dimension "without proof by 
Bray and Blythe Q to calculate an upper bound for 
the survival probability for d < 2, although it was al- 
ready proved in 10] in the context of incoherent exciton 
quenching. More recently, the principle was proved by 
Bray, Majumdar and Blythe Almost simultaneously 
but slightly earlier, Moreau et al. ^2 proved the Pascal 
principle for a rather general class of random walks on d- 
dimensional lattices. These latter authors are responsible 
for the appealing name now bestowed on the principle. 

The purpose of this paper is to extend the procedure 
and results of Bray and Blythe |^ , which are valid for a 
Brownian diffusive particle and Brownian diffusive traps, 
to situations in which the particle and traps move sub- 
diffusively. Anomalous diffusion of a particle is usually 
characterized by its mean squared displacement x{t) for 
large t: 



{x\t)) 



2K-, 



r(i + 7) 



(1) 



Here is the (generalized) diffusion constant and 7 
is the anomalous diffusion exponent. Ordinary Brow- 
nian diffusion (7=1, Ki = D) follows Fick's second 
law, (^x^{t)^ (X t. The process is called sudiffusive when 
< 7 < 1. Subiffusive processes are ubiquitous in na- 
ture [Hill El E milil, and are particularly useful 
for understanding transport in complex systems 0, Il9| . 
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The problem considered in this paper is a special case 
of a broad class of reaction- subdiffusion processes that 
have been studied over the past decades. One approach 
that has been used to study these processes is based on 
the continuous time random walk (CTRW) with waiting- 
time distributions between steps that have broad long- 
time tails and consequently infinite moments, •ip{t) ^ 
t-i-T for t ^ oo with < 2< 1. Using the CTRW 
formalism, Blumen et al. |23, |21| considered a variety of 
reactions including the trapping problem A+B{steLtic) — > 
i3(static), the target problem A(static) + B ^ B, and 
the bimolecular reactions A + A ^ (d and A + B 0. 
The moving particles were modeled as continuous-time 
random walkers with long-tailed waiting-time densities. 
Quite recently, Sung and Silbey have used the CTRW 
model to study the dynamics of particles that react at a 
boundary. A CTRW approach has also been applied by 
Seki et al. l23!| to study the kinetics of the recombination 
reaction in subdiffusive media. 

Another approach is based on the fractional diffusion 
equation, which describes the evolution of the probability 
density P(x, t) of finding the particle at position x at time 
t by means of the fractional partial differential equation 
(in one dimension) pll23| 



di 



P{x, t) = oA 



92 



(2) 



where is the generalized diffusion coefficient that ap- 
pears in Eq. J^l, and qD\~^ is the Riemann-Liouville 
operator, 

1 f) 



dr-^^. (3) 



r(7) at Jo {t - t) 

Sung et al. directly addressed this problem with a 
fractional diffusion equation approach. Seki et al. [2^ 
went beyond the CTRW model and derived a fractional 
reaction-diffusion equation for the geminate recombina- 
tion problem, and one finds some disagreement between 
the assumptions and results in this work and that in |25j| . 
The fractional diffusion approach has recently been used 
to get exact solutions for two types of one-dimensional 
trapping problems: the so called one-sided problem, in 
which all the traps lie on one side of the particle, and 
the two-sided problem, in which the traps are located 
on both sides of the particle (this is the traditional or 
standard version of the trapping problem) |26i |. The 
fractional-diffusion approach has also been employed to 
study other bimolecular reactions between subdiffusive 
particles. In particular, the annihilation A + A $ and 
coa gula tion A-t- A — > Aoi subdiffusive particles was stud- 
ied [231 by means of a fractional generalization of the 
interparticle distribution function method 4]. The evo- 
lution of reaction-subdiffusion fronts ior A + B C re- 
actions, where both A and B move subdiffusively, is also 
amenable to analysis by means of the fractional diffusion 
approach 28]. Other recent work on fractional diffusion 
and CTRW models of subdiffusive reacting particles can 
be found in a number of references 1291. 



In this paper we implement the fractional diffusion 
equation approach to study the one-dimensional trap- 
ping problem in the long-time regime for subdiffusive 
(or diffusive) particles that move among a distribution 
of non-static traps. The traps can be either subdiffusive 
or diffusive. For this purpose, we generalize the ideas of 
Bray and Blythe ^] . Recent contributions to the A-{- B 
problem based on the fractional diffusion equation ap- 
proach [2^ |26| share the simplifying characteristic that 
the reaction takes place between a static particle (or fixed 
boundary) and a subdiffusive particle. The present work 
differs from those in that all the reacting particles (in- 
cluding traps) are (sub)diffusive and, moreover, the dif- 
fusion constant and the anomalous diffusion exponent of 
each species may be different. Some of this work has been 
presented in 
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In some cases, asymptotic anomalous diffusion behav- 
ior can be found from corresponding results for normal 
diffusion with the simple replacement of t by t''' . This can 
be understood from a CTRW perpective because the av- 
erage number of jumps n made by a subdiffusive walker 
up to time t scales as (n) ~ , and in many instances the 
number of jumps is the relevant factor that explains the 
behavior of the system. The simple replacement result is 
evidence of "subordination" (see Sees. 5 and 7.2 of ^21] ) . 
However, there are other instances where the behavior 
of subdiffusive systems cannot be found in this way. A 
simple example is the survival probability of subdiffusive 
particles in the trapping problem (see Sec. 5 of [IJ). 
In particular, for systems where each species has a dif- 
ferent anomalous diffusion exponent, such a replacement 
becomes ambiguous. This is the case for the problem 
considered here. 



Bray and Blythe obtained the asymptotic survival 
probability of a diffusing particle in a sea of diffusing 
traps by calculating an upper and lower bound that con- 
verge asymptotically. We follow their procedure for sub- 
diffusive particle and traps, with partial success. While it 
is possible to obtain convergent bounds in most anoma- 
lous diffusion exponent regimes, this procedure does not 
work in all regimes. In particular, the bounding proce- 
dure encounters difhculties when the particle A diffuses 
normally and the traps are "too slow" but not static. 



Our paper is organized as follows. In Sec. ^ we extend 
the proof of the Pascal principle in one dimension to the 
case in which both the particle and the traps move subd- 
iffusively and calculate the upper bound for the survival 
probability. In Sec. lIIII we present an alternative calcula- 
tion of the upper bound, which is particularly helpful for 
the calculation of the lower bound in Sec. IIVI The sur- 
vival probability is established, when possible, in Sec.lVl 
Section IVII presents a compendiary of results and some 
comments on open problems. 
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II. THE PASCAL PRINCIPLE 

The "Pascal principle" of random walks says that the 
best survival strategy for a random walker A surrounded 
by a random sea of trapping random walkers B is to stand 
still. Here we extend the proof of the Pascal principle in 
one dimension to the case in which both the particle and 
the traps perform subdiffusive random walks. 

As did Bray et al. we consider a finite volume V 
containing N = pV traps B initially distributed at ran- 
dom, and a single A particle initially at the origin. The 
trajectory of the A particle is z{t). Bray et al. write 
the survival probability of A as P{t) = exp{— /Lt[z(t)]}, 
where the trajectory-dependent functional jj. is to be de- 
termined. To find this functional they derive an equation 
for it by calculating, in two ways, the probability density 
to find a B particle at the point z(t) at time t: 



dt' ii{t')G[z{t),t\z{t'),t']+ PP W)) ■ (4) 



That the left side is this probability density is obvious. 
On the right side ^i{t')dt' = {-P/P)dt' is the probabil- 
ity that a B particle intersected A in the time interval 
{t\ t' + dt') for the first time, and the propagator G is the 
probability density for this particular B to be at z(t) at 
time t. We have slightly augmented Bray et al.'s funda- 
mental starting point by including the second term on the 
right hand side, which is the probability that the first in- 
tersection of a -B with A occurs at time t and not before. 
This term is ultimately unimportant because it decays 
much more quickly than the first term and so we omit 
it henceforth, but it is satisfying that the fundamental 
equation now holds at all times. 

Proof of the Pascal principle requires us to show that 
the trajectory z{t) — for all t gives the smallest pos- 
sible value of iJ-[z]. For this trajectory the fundamental 
equation is 



p= / dt'fio{t')G[zit)^0,t\z{t')^0,t'], (5) 
Jo 

where /io(i) = fj,[z{t — 0)]. The propagator 
G {z{t) = 0, t\z{t') = 0) is a function only oi t — t', which 
we abbreviate as G{t — t'), so Eq. ||SJ) is a convolution. 
Denoting the Laplace transform of f{t) as /(s) yields for 
the transform of ||SJ) 



p{s) = s/io(s)G(s) 



(6) 



where we have noted that pio{t' = 0) = because the 
initial survival probability is unity. It follows that 



Ao(s) 



sG{s)' 



(7) 



Inversion for the subdiffusive case in one dimension gives 

t-7/2 



Git) 



H 



1,2 



(1-7/2,7/2) 
(0,1), (1/2,1) 



y47?^r(i-i)' 



(8) 



where 7 is the anomalous diffusion exponent for the traps, 
K~f is the associated generalized diffusion constant, and 
is Fox's iJ-function [siIIb^ . whose value at the given 
arguments we have used to write the last equality. The 
Laplace transform of G{t) is 



G{s) = 



The inverse of Eq. Q then immediately follows: 



r(i + i)' 



(9) 



(10) 



To prove that z{t) — gives the global minimum of 
/i[z(t)] we again follow Bray et al. and write /i = //q + Mi 
in Eq. Q (without the last term): 



p= / dt' [fi^it') + fii{t')]G[z{t),t\z{t'),t'] 



(11) 



Adding and subtracting {t - t')-'^/'^ / ^AK^T (l - ^) to 
G in the integrand allows cancellation of the left side of 
Eq. Hll() against one of the terms on the right, leaving 



= - / dt' 



R 



(1 - (t - 1'^/^ 

1 Ai(i') 



dt' 



where 

R = 



r(i-i).^2.o 



1,2 



(1-7/2,7/2) 

(0,1), (1/2,1) 



(12) 



(13) 



An explicit expression for /xi is obtained by Laplace 
transforming this result, solving for p.{s), and inverting. 



Mi(<) = 



dt^ 



r(i)r(i-i)7o {t-t,Y 

/*i Rfijt^) 
X / dt2 -, / ,^ . 

h ih-t^r^^ 



(14) 



The rest of the argument follows exactly as in Bray 
et al. 0, [33 ■ Since R > and /i > (because 
P{t) is a nonincreasing function of time), it follows 
that fJ-i{t) > for all paths z(t), with equality when 
z{t) = for all t. The survival probability of our 
A particle averaged over all possible trajectories z is 
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P{t) = (e-''(*));, = (e-''o(*)-''i(*))^, and the above proof 
shows that {e->^oW-t^iW) ^ < e-A'o(*). Therefore, 



Pu{t) = exp [-Mo(0] = exp 



r(i + i) 



(15) 



is a strict upper bound for the survival probabihty P(i), 
thus extending the proof of the Pascal principle to subd- 
iffusive particles in one dimension. 



III. ALTERNATIVE CALCULATION OF UPPER 
BOUND FOR THE SURVIVAL PROBABILITY 

The survival probability Pu{t) of a static particle A 
surrounded by a distribution of mobile traps (the "target 
probler n" ) was considered as early as 1986 using a CTRW 
model |2l|, and more recently for the three-dimensional 
case by means of a fractional diffusion approach |^. In 
the previous section we obtained the result ((TH|l following 
the approach of Bray et al. ^Jl- Here we re-calculate 
Pu {t) bygeneralizing the original approach of Bray and 
Blythe Q to the subdiffusive case because it provides 
results useful for the calculation of the lower bound in 
the next section. 

Consider a target A of size 2C centered at the origin, 
and let Qi{t\y) be the probability that the trap initially 
placed at 2/ > £ has not reached the end of the target 
at y = £ by time t. Then, in terms of the Fox's H 
function ISJ. 



Qi(t|y)-l-i/i\" 



= 1-H 





(1,7/2) " 




(0,1) 



(16) 



For 7 — > 1 the Fox's H function becomes the comple- 
mentary error function (with Ki = D), and the ordinary 
Brownian motion result is recovered, 



Qi{t\y) = l-erfc 



y~c 



1 



(17) 



Next consider iV independently diffusing traps that are 
initially placed at random in the interval C < y < L + 
R. Here and henceforth 2R is the size of the system, 
which we will take to infinity at appropriate points in the 
calculation. The probability QN{t) that the stationary 
target A has survived up to time t is 



A' 



dy, ll-H 



yi- 1^ 



1 

'r 



C+R 



dyH 



N 



(18) 



or, in terms of the density p = N/R of traps, 

rR 



Qoo(i) = lim ■ 



exp 



-I \ pR 



y 



dzH\ 



(19) 



Note that the result is independent of the size of the 
target. We need to evaluate the integral 



dz Hl° 



(1,7/2) 
(0,1) 



(20) 



which can be done from the properties of the Fox's H 
function [^J- One finds that 



rrlO 

^11 



(1,7/2) 
(0,1) 



trio 



(1+7/2,7/2) 
(0,1) 



But 



and 



^11 



oo 



(1-^7/2,7/2) 
(0,1) 



Trio 



(1 + 7/2,7/2) 
(0,1) 



1 



r(i + 7/2)' 



so that I-y = l/r(l + 7/2). Therefore, 



Qoo(i) = exp 



"r(i + 7/2) 



(21) 
(22) 

(23) 
(24) 



This is the survival probability of the target when all the 
traps are located to its right. When the traps are located 
on both sides of the target, the survival probability of the 
target is the square of Eq. (|24|l : 



Pu{t) - Qloit) - exp 



"r(i + 7/2) 



(25) 



This result, which is of course identical to Eq. (|15|l . is the 

upper bound on the survival probability of the moving 
particle. 

Incidentally, as is well known, the survival probability 
for the target problem is related to the distinct number 
of sites S{t) visited by a trap up to time t ^21. .34] . 



Pu{t) 



^ e-p{5(t)>+- 



(26) 



where the dots represent higher moments that decay 
more rapidly with time. Comparing this expression with 
Eq. (|25|l . one finds that the asymptotic average value 
{S{t)) of the territory explored up to time i by a sub- 
diffusive walker with generalized diffusion coefficient Kj 
and anomalous diffusion exponent 7 is 



{S{t)) 



(27) 



r(i + 7/2)- 

This result agrees with that found by Yuste and 
Acedo 01 using a different approach. 
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IV. LOWER BOUND FOR THE SURVIVAL 
PROBABILITY 



the subdifFusive case, we note that for large arguments 
{z ^ 1) the Mittag-Leffler function has the expansion 



Let PL{t) be the probabihty that the mobile parti- 
cle A remains inside a box of size C and that all the 
traps remain outside this box until time t (we distin- 
guish between C, the size of the box, and L, which de- 
notes "lower bound"). When this happens, the particle 
A survives. It is clear that PL{t) is a lower bound for 
the survival probability P{t) of interest because there 
exist many other trajectories involving the simultaneous 
presence of the particle and traps within the box C that 
allow the particle A to survive. This lower bound was 
first calculated for diffusive particles and traps by Red- 
ner and Kang |35| and further considered (and in some 
cases corrected) in 0, |3^ . 

The probability PL{t) is itself the product of three 
probabilities: 

1. The probability Qi that at t = the box of size £ 
contains no traps: 



Qi 



(28) 



2. The probability Q2 that no traps enter the box of 
size C up to time t: 



Qi = exp 



r(i + 7/2) 



Pu{t). (29) 



Note that it is the derivation of the previous section, 
which explicitly shows this probability to be independent 
of the size of the box, that allows us to write this result. 

3. The probability that the particle has not left the 
box of size C up to time t. We proceed to evaluate this 
quantity. 

Let VF(a;,t) be the probability of finding the particle 
A at position x at time t if it was at position x = 
at time t = Q and there are absorbing boundaries at 
X = ±£/2. Solving the fractional diffusion equation by 
means of separation of variables ,13] one finds 

m^,t)4f:(-l)"sini^!i±lM^i±^ 

ri=0 

X E^, (-if;,(2n + Ifir'^t^'/L'^ , (30) 

where K'^, and 7' are the generalized diffusion constant 
and the anomalous diffusion exponent of the particle A. 
Therefore, 



(-1) 



m+l 



SO that 



™=i r(l - 7'm) 



r(i — 7'm) 



(32) 



(-1)^ 



7 



Therefore, for t ^ 00 one finds 
1 £2 



O 



£2 



(33) 



(34) 



Consequently, a lower bound on the survival probability 
of the particle A is 

PL{t) ^QiQiQs 



-e exp 



L r(i + 7/2) 

£2 



8r(i - 7') K'tr 



1 + 



(35) 



It can easily be ascertained that this expression is max- 
imal when C = C* = 2/ p (independent of time), i.e., 
Pdt) < PL-it) with 



PL-{t) 




(36) 



This then is our best lower bound for the survival prob- 
ability P{t) of a subdiffusive particle. 

When the particle A diffuses normally, Eq. H31|l be- 
comes 



TT 2n + 1 

n=0 



with D' = K[. For long times 



C/2 



-C/2 



W{x,t)dx 



i-iy 



A 00 

-y 

TT 2n 

n=0 



1 



-if;,(2n-H)27rV/£2 . (31) 



Next we distinguish two cases in the handling of the 
sum in Eq. H31|l : first we deal with a subdiffusive particle, 
and subsequently with an ordinary diffusive particle. In 



Q3 exp [-^D'nh/C^] , i > 1 



so that 



Pdt) =QiQ2Qz 



^ -oC 

j—e exp 

TT 



r(i + 7/2) 

exp [-D'TTh/L^] 



p^K-yt-y 



(38) 



(39) 



6 



for t ^ 1. This lower bound can again be maximized 
by optimizing the value of £. The optimal value is ^ 

1/3 

ime dependent), so that 



PL{t) < PL'{t) 

4 

= — cxp 



"r(i + 7/2) 



1/3 



(40) 



Note that the dominant term inside the bracket depends 
on the value of 7, the anomalous diffusion exponent for 
the traps. We distinguish three cases: 

1. Traps with 2/3<7<l. In this case, for t ^ 1, 



r(i + 7/2) 



Nl/3 



SO that 



PL'{t) = -exp 



r(i + 7/2) 



(41) 



(42) 



2. Traps with 7 = 2/3. Now 
PL{t) < PL'{t) 



4 

— exp 



r(4/3) 



(43) 



that is, the second contribution in the exponent in 
Eq. (|40|l is of the same order as the first and must thus 
be retained. 

3. Traps with < 7 < 2/3. Now the second term in 
the exponent of Eq. (|4U|I is dominant: 



Pi.(f) = -exp -3(7rV^£''V4) 



Nl/3 



(44) 



In the next section we examine our upper and lower 
bound results to establish the behavior of the survival 
probability of A whenever possible. 



V. SURVIVAL PROBABILITY 

We now combine our upper and lower bound results. 
Recall that the label and exponent 7 is associated with 
the traps and 7' is associated with the particle A. The 
upper bound on the survival probability is in all cases 
given in Eq. (|25|l . but the lower bound depends on the 
anomalous diffusion exponent of the particle. We distin- 
guish the following cases: 

1. Subdijfusive particle (0 < 7' < 1) and diffusive or 
subdiffusive traps (0 < 7 < 1). The lower bound is given 



in Eq. so that PL-it) < P{t) < Pu{t) leads to 



< - 



In P(t) 

r(i + 7/2)- ^/^W^ 
2 



< 



r(i + 7/2) 

2 In 



p^K'ti' +2 + ln[2r(l-7')] 



^/p^K^ 



O 



p^^,ty 



(45) 



p'^K'tt' < ^ p^K^t"^ and the upper 



For t — > 00, In 

and lower bounds converge asymptotically. We therefore 
arrive at the explicit asymptotic survival probability 



P(t) - exp 



r(i + 7/2) 



p^K^t-^ 



(46) 



for < 7 < 1 and < 7' < 1. Note that for 7 = 1 we 
recover the normal diffusive result obtained earlier j^. 
A noteworthy result here is that the survival probability 
depends only on the exponent 7 that characterizes the 
traps and not on 7' that characterizes the particle. This 
is interesting vis a vis the subordination issue. 

2. Diffusive particle (7' = 1) and subdiffusive traps 
with 2/3 < 7 < 1. The bounds here are 



r(i + 7/2) 



< 



< 



InP(i) 
2 

r(i + 7/2) " 



_ > 2/3 ^,1/3 



2p 



K. 



1/2 

7 



(47) 



and the asymptotic survival probability is again given by 
Eq. ggi- 

3. Diffusive particle (7' = 1) and subdiffusive traps 
with 7 = 2/3 (marginal case). Now PL-it) < P{t) < 
Pu{i) leads to the more ambiguous inequalities 



r(4/3) 



< 



< 



InP(t) 
2 



r(4/3) 



^^2/3 ^,1/3 

2p, 



Ill ■ 



(48) 
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The bounding procedure is therefore not able to predict 
the value of the prefactor A in P(t) = exp(— At^/'^), but 
the asymptotic behavior — In P(t) cx i^^'^ is evident. 

4. Diffusive particle (7' — 1) and subdijfusive traps 
with < 7 < 2/3. The bounds here are also given by 
Eq. (|47|l . so that the bounding procedure is not able 
to determine the asymptotic behavior of P{t) at all for 
this case. We are not even able to assert the asymptotic 
stretched exponential form P{t) exp {-Xt^). 
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VI. PANORAMA AND DISCUSSION 

Bray and Blythe Q have calculated the asymptotic 
survival probability of a diffusive particle A in a randomly 
distributed sea of diffusive traps B in one dimension, and 
have determined the precise value of the coefficient A in 
the classic result P{t) ~ exp (— At^/^) first obtained by 
Bramson and Lebowitz |a,l3|- Within some constraints, 
we have generalized this result to the case where one or 
both of the species move subdiffusively. Our particle A 
is characterized by the anomalous diffusion exponent 7' 
and generalized diffusion coefficient Ky , and the traps by 
7 and K^. In the process of this generalization, we have 
extended the proof of the Pascal principle, that the best 
survival strategy of a particle in a sea of moving traps is 
to remain stationary, to the case of particles and/or traps 
that move subdiffusively. These may be the first results 
in the literature involving two subdiffusive species with 
different anomalous diffusion exponents. 

When both species are subdiffusive (7 and 7' both 
smaller than unity), the survival probability is indepen- 
dent of 7' and determined entirely by the subdiffusive 
properties of the traps, cf. Eq. When the par- 

ticle moves diffusively (7' = 1), on the other hand, 
we are unable to unequivocally determine the coeffi- 
cient A for all cases using this procedure. If the traps 
move sufficiently rapidly (2/3 < 7 < 1) then the result 
Eq. H46|) is still valid. Note that this reduces to the Bray 
and Blythe result when 7 = 1. The case 7 = 2/3 is 
marginal in the sense that we can establish the behavior 
P{t) ^ exp (— At^/-^), but are not able to determine the 
constant A. Note that this particular time dependence 
of the survival probability is the same as the classic re- 
sult for the survival probability of a diffusive particle in 
a sea of immobile traps If the traps are too slow 

("strongly subdiffusive"), < 7 < 2/3, we are not able 
to determine even the time dependence of the survival 
probability on the basis of this approach. However, since 
we find the same stretched exponential behavior when 
7 = 2/3 and when 7 = 0, a conjecture as to the behavior 
throughout this slow trap regime might be appropriate. 
The conjecture is that the survival probability decays as 
P{t) ~ exp(-Ati/^) in the entire regime < 7 < 2/3. 

We thus find that in so far as one can even think of 
some sort of subordination principle (and whether such 
thinking is appropriate here is debatable), it is deter- 
mined by the behavior of the traps, i.e., by the replace- 
ment of t by f. Even in the range of exponents where 
this is possible, it is only possible for the main asympotic 
contribution to P{t) but not for the correction terms to 
the leading asymptotic term. 

It is interesting to note that the value £* — 2/p that 
maximizes the lower bound of the survival probability 
for a subdiffusive particle A does not grow with time. 
This implies that finite particle size effects could become 



relevant with increasing density p. This is completely 
different from the case of a Brownian particle A, since 
the growth C* cx t^^^ now suppresses such finite size con- 
tributions for any given density. 

At this point we inject a digression that is relevant not 
only to our analysis but also to the original work of Bray 
and Blythe They assumed that the particle A is ini- 
tially surrounded by a random (Poisson) distribution of 
mobile traps, an assumption also made in our explicit 
analysis, cf. Eq. (|28|l . On the other hand, if at the start 
of the observations {t = 0) the process has already been 
taking place for some time — r (i.e., if the process started 
at some time r in the past), then it is known that the 
distribution around the surviving particles at time t = 
is not of Poisson form. Those particles that initially had 
close-by traps are more likely to have been trapped al- 
ready than those that did not, so that those particles 
that have survived are surrounded by a region of fewer 
than average traps (sometimes referred to as a "gap"). 
Bramson and Lebowitz arrive at the conclusion that the 
configuration of B particles is nevertheless dominated by 
a Poisson random measure 0. In [lOl we confirm that 
for any finite r the gap does not affect the asymptotic 
survival probability results of Bray and Blythe. The de- 
tailed nature of the gap is different in the diffusive and 
subdiffusive cases, and unknown in the latter. However, 
we conjecture that it is no more pronouced in the subd- 
iffusive than in the diffusive system, and that it does not 
affect our results cither. 

Our own results of course leave a number of questions 
unanswered. One obvious question concerns the marginal 
role of the trap exponent 7 = 2/3 when the particle is dif- 
fusive. Why is this a marginal exponent? A connection 
between this critical value and the fact that for a Brow- 
nian particle the length that maximizes the lower bound 
of the survival probability grows as C* ^ t^^^ seems plau- 
sible, but the conceptual basis for such a relation is not 
clear. 

The most pressing and intriguing puzzle to resolve is 
that of calculating the survival probability when the par- 
ticle A is diffusive (7' = 1) and the traps are strongly 
subdiffusive (0 < 7 < 2/3). While we conjecture that 
the survival probability in this regime decays as P{t) ^ 
exp(-Aii/3), the upper and lower bounds in this case do 
not have the same asymptotic time dependence so we 
are not able to test these conjectures on the basis of the 
procedures used in this paper. 
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